Intercommutation of Semilocal Strings and Skyrmions 

Pablo Laguna 

Department of Astronomy & Astrophysics, IGPG, CGWP 

Penn State University, University Park, PA 16802 

Vishnu Natchu and Richard A. Matzncr 

Center for Relativity, Department of Physics 

The University of Texas at Austin, Austin, TX 78712 



O 
O 
(N 

< 

oo 
> 

o 
\o 
o 

-f— > 
I 

a*: 



X 



Tanmay Vachaspati 

CERCA, Department of Physics, Case Western Reserve University, 

10900 Euclid Avenue, Cleveland, OH 44106-7079 

We study the intercommuting of semilocal strings and Skyrmions, for a wide range of internal 
parameters, velocities and intersection angles by numerically evolving the equations of motion. We 
find that the collisions of strings and strings, strings and Skyrmions, and Skyrmions and Skyrmions, 
all lead to intercommuting for a wide range of parameters. Even the collisions of unstable Skyrmions 
and strings leads to intercommuting, demonstrating that the phenomenon of intercommuting is 
very robust, extending to dissimilar field configurations that are not stationary solutions. Even 
more remarkably, at least for the semilocal U(2) formulation considered here, all intercommutations 
trigger a reversion to U(l) Nielsen-Olesen strings. 



Physical systems can have a large number of sponta- 
neously broken symmetries, some of which may be local 
(gauged) while others may be global. In cases where 
the symmetries are only partially gauged, or gauged 
with unequal strengths, semilocal strings may be present 
UlSSul- If the gauge couplings are widely disparate, 
the strings may also be stable. Semilocal string solu- 
tions exist in the standard electroweak model, but the 
SU(2) and U(l) gauge couplings are close enough that 
the solutions are unstable 0, Q- Recently twisted and 
current-carrying semilocal string solutions have been dis- 
covered that may be stable for other values of parameters 
0,@]- Semilocal strings arise in supersymmetric QCD if 
the number of flavors is larger than the order of the sym- 
metry group [3- In string theory, they can arise as field 
theoretic realizations of cosmic D-strings 10], in which 
case they may play a role in cosmology. The formation of 
semilocal strings in a phase transition has been studied in 
Refs. [Ulli^l) in a cosmological setting in Refs. |l3lll4j. 
and in superstring cosmology |15j . 

In the case of topological gauged [T|| E3, 13 an d global 
[l9l| U(l) strings, numerical evolution of the field the- 
ory equations show that colliding strings intercommute 
(Fig. for almost [2(j any scattering angle and veloc- 
ity. Thus the probability of intercommutin g is unity for 
topological strings. Qualitative arguments [n| HH have 
been given in an effort to understand intercommutation 
of topological strings. In the context of string theory cos- 
mic strings and QCD strings, intercommuting probabili- 
ties have been calculated in Refs. |23l l24| within various 
approximations. 

In this paper we study the interactions of semilocal 
strings, which are topological structures embedded in 
the solution space of the partially gauged theory, and 
the related semilocal Skyrmions which are fat and fuzzy 
objects (still with one "long", i.e. stringlike dimension) 



on which strings can terminate, much like strings ending 
on monopoles. (We shall show below (Figs ^ how- 
ever, that the termination is dynamic, moving along the 
string.) The scalar field in a Skyrmion can be arbitrar- 
ily close to its vacuum expectation value everywhere, in 
contrast to a string in which the scalar field magnitude 
vanishes (which is a nonvacuum value) on some curve. 
We review below that the embedded topological strings 
have pre cisely the dynamics of U(l) (Nielsen-Olesen 
(NO) 22]) strings. Thus semilocal strings will intercom- 
mute on collision with each other. However, the inter- 
action of different objects such as semilocal strings with 
Skyrmions, or Skyrmions with Skyrmions, is a new con- 
sideration. Remarkably we find that intercommuting is 
extremely robust - strings intercommute with Skyrmions, 
and Skyrmions intercommute with other Skyrmions. In- 
deed, even in the parameter regime where the solutions 
are unstable, we see intercommuting. This establishes, 
for the first time, that intercommuting is not particular 
to string solutions, but can happen more generally be- 
tween any two field configurations. Even more remark- 
ably, at least for the semilocal U(2) formulation consid- 
ered here, all intercommutations trigger a reversion to 
U(l) NO strings. 

The semilocal model we shall study is the bosonic part 
of the electroweak model with weak mixing angle W — 
7r/2. The Lagrangian is: 
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L = |(0„ - iY u )<S>\ 2 - -Y^Y^ - £($t$ - 1 



(1) 



where $ T = (<p,(x) is an SU(2) doublet, Y~ M is a U(l) 
gauge potential and Y^ = d^Y v — d^Y^ its field strength. 
The only parameter in the model is j3. 

The semilocal model has [51/(2) x U(1)y]/Z 2 symme- 
try and only the U(1)y is gauged. If $ T gets a vacuum 



expectation value along the direction (0, 1) the symme- 
try is broken down to U(1)q where Q — (T 3 + Y)/2 = 
diag(l, 0), T 3 = diag(l, -l)/2, and Y = diag(l, l)/2. 

The NO string is an embedded string solution to 
Eq.(UJ): Hli$ T = (0,fe lf> ) and Yg = (v/r)eg in 
cylindrical coordinates (r, 9, z) with f(r) and w(r] profile 
functions identical to those for the NO string 22]. Fur- 
ther, the dynamics and interactions of fields living in the 
embedded subspace - lower component of $ and Y^ - are 
given exactly by the Abelian-Higgs model and so these 
behave in exactly the same way that NO strings of the 
Abelian-Higgs model do, and intercommute regardless of 
collisional velocity and angle |la . fl7J fl8| . However, in this 
subspace of the field theory, excitations that can cause 
the string to terminate in Skyrmions are absent. So to 
study the consequences of semilocality for string interac- 
tions, we need to go beyond the embedded Abelian-Higgs 
model. 

It has been shown that the NO string, as embedded in 
this model, is stable for (3 < 1, unstable for j3 > 1 and 
neutrally stable for j3 = 1 ■ This means that if we start 
with a string-like configuration that involves degrees of 
freedom outside the embedded Abelian-Higgs model, it 
will relax into the string solution for [3 < \ and go further 
out of the embedding space for (3 > 1. Hence the most 
suitable choice of parameter for studying the effects of 
semilocality on the interactions of strings is precisely at 
[3 = 1 and this is the value that we adopt for most of our 
analysis. 

At (3 = 1 there is a class of solutions that are degener- 
ate in energy and that can be written as |27| : 
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where u = In |$| 2 is the solution to: 

V 2 w + 2(l-e u ) = V 2 ln(r 2 



u(r;\q \) 



W\ 
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and u — > as r — > oo. The complex parameter qq is a 
coordinate on the space of solutions. At qo = 0, the solu- 
tion is the embedded NO string. For other values of \qo\ 
the magnetic field and the energy density is more spread 
out and the solutions are called "semilocal Skyrmions" 
|13| or, generically, semilocal strings with particular val- 
ues of qo. The phase of qo labels a Goldstone mode - for 
example, the phase can be pulled out as an overall fac- 
tor of the field <i>, and be compensated for in the lower 
component by shifting 8. The relative phases of qo for 
two strings (called the relative "color" in [23) can still 
be significant but the phase difference is expected to be 
quickly radiated away as the two strings collide |2g, so 
it is sufficient to take go to be real. In fact, as described 
below, our present simulations show that qo — > in the 
aftermath of the collision. 

We now set up two semilocal strings with (real) qo 
values q^ and <? ■ The strings are then boosted so that 
they approach each other at relative velocity, V (in units 
of c), and angle O. Numerical evolution of the equations 



then tells us whether the strings intercommute or pass 
through. 

Initial data and boundary conditions are based on suit- 
able superpositions of single, boosted, static string solu- 
tions. Static solutions $ T = (g, fe ) and Yg = (v/r)e$ 
for semilocal strings with (3=1 can be found from: 
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where /, g and v are functions of r only and g — q f/r. 
To set up initial data and boundary conditions one could 
numerically solve Eqs. Q and JSJ and construct a table 
to be used for the superposition of string solutions. We 
choose instead parametrized approximate solutions to the 
equations. That is, 
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where f NO = [tanh (r Rf/2)] B f ,R f = l + (a 2 -l) e-' r / r ° 
and Ef = a 3 + (1 — a4)e~ r > r ° with similar expressions 
for v. The parameters in these expressions are found by 
minimizing the total energy for a single string. 

Note that the fields for semilocal strings with non-zero 
go approach their vacuum values only as power laws (r~ 2 
for $, and r -4 for the magnetic field strength). This is 
to be contrasted with the exponential approach for NO 
strings (which is the qo = case). Nonetheless these 
polynomial fall-offs are still fast enough that we antici- 
pate being able to construct initial data by superposing 
fields for widely separated strings, with minimal change 
in the total energy. The scheme is remarkably straight- 
forward. (An early reference to a similar scheme for NO 
strings is Abrikosov |29|. See also Matzner J16J.1 

Next we consider how to superpose the scalar fields of 
two well-separated strings, labeled by i = 1,2. Let the 
two string scalar field configurations be $f = (4>i,<7i) 
and $>2 — (02, 02), respectively. Since the <pi vanish as 
a power law at infinity, and Oi go to unity, we superpose 
the scalar field string configurations using the scheme 
<1> T = (0i + 02 ,0-1 02). Because the axial gauge field Y^ 
satisfies an equation which is linear, the superposition of 



Y„ is taken to be Y„ = Y, ' 



Y L 



(2) 



M 10 »av,ii u" «v> ijB - j m -r ± fi . 

The actual data are set as follows. With the energy 
minimization procedure described above, we start by con- 
structing two single string solutions with parameters qo 
and q^ in coordinates (x,y, z). The strings are parallel 
to the z axis and pass through the origin. We hence- 
forth assume that Y^ is written in the Lorentz gauge: 
d^Y^ = 0. This gauge is linear and Lorentz invariant. 

(2) 
The string q is then rotated around the x-axis by an 

angle 6, with O = (parallel strings) and O = 7r (an- 

tiparallel strings). Next, both strings are off-set along the 



x-axis a distance x 



(i) 



f2l 

D/2 and x = D/2, respec- 
tively. Finally, the strings are Lorentz boosted toward 





(a) Initial configuration 



(a) Initial configuration 





(b) Late configuration 



(b) Late configuration 



FIG. 1: Field a for a collision of go = 1 ("top") and go = 3 
scmilocal strings, at Q = 90°, with collision velocity V = 0.9. 
The contours arc 1 — |<r| 2 at values 0.1,0.2,0.5. Note that 
l°"| 2 + \4>\ 2 — 1 is the vacuum. The strings reconnect. 



each other with velocity V, so that the string intersec- 
tion will occur in the middle of the computational cube. 
Initial time derivatives are computed to second order in 
the timestep by repeating the data construction process 
at time At/2 into the future and At/2 into the past and 
taking 1/Ai times the difference of the two field configu- 
rations. 

The boundary data are naturally set as part of the ini- 
tial setup process. In subsequent time steps, the bound- 
ary values are computed as if this procedure were re- 
peated, with the strings advanced at the boundary as 
if in free motion with the initial velocity. This means 
that even after the collision in the center of the cube, the 
boundaries continue to act as if the strings had passed 
through one another, regardless whether they did in fact 
pass through one another, or reconnected. However our 
computational domains are large enough that the region 
where the strings collide is causally disconnected from 
the boundaries of the computational domain. 

Our numerical code solves the field equations of motion 



FIG. 2: Coincident configurations of the field <f> in the collision 
of FIG 1. The contours are |0| 2 at 0.2,0.3,0.5. The value 
<j> — is an NO string. It can be seen that the configurations 
reverts to NO form after the collision. This is in spite of 
the fact our boundary conditions hold the boundary values 
of the strings in their original go = 1 and go = 3 state. The 
transition from the initial value of <j> to the value of zero means 
we have two NO strings near the center terminating on the 
(outward moving) ends of the g 7^ semilocal strings. 



for $ T = (<f>, a) and Y^ as a coupled set of first-order 
differential equations in time: 



dtVA = VtA-Y'tA-d^U 

-2e AB (y t r ?B + y J 9^ B ) 
8 t K A = V 2 a A -Y 2 a A -d aA U 

-2e AB {Y t K B +Y l d l n B ) 
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where d t <f) A = r? a, d t a A = n A and <9 t Y" M = W^. The 
subscripts A, B = 1, 2 denote respectively real and imag- 
inary parts, e A B is the completely antisymmetric tensor 
in two indices with ei2 = 1 and U — [3((t> 2 + a 2 — l) 2 /2. 
We use a second order accurate, both in space and time, 



discretization. The temporal updating is done via a stan- 
dard leap-frog method. 

Our first set of numerical experiments were done with 
strings of the same internal parameter go ( % = % )• 
We first verified that go = strings (NO strings) 
reconnect, as had been found previously. We also 
considered collisions between q = 1 strings. Just 
like NO strings, these qo = 1 semilocal strings al- 
ways intercommuted for the set of velocities V G 
{0.1,0.5,0.9} = {V} and all intersection angles tested: 
{30°, 45°, 60°, 90°, 120°, 150°} = {©}. Next, we consider 
Our simulations now depend on 



,1 (2) / (1) 

the case q f= q 
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four parameters: V, O, q$ ' and Qq — Qq ■ To reduce 
the parameter space, we only considered g^ = 1, and 

q ( 2) G {0.0,0.5,0.9,1.0,1.1,2.0,3.0}. In these cases also 
we find that for all velocities V G {V} and all angles 
O G {©}, reconnection always occurs. The dynamics 
involved in the reconnection are interesting. In particu- 
lar, the value of go for both strings appears to relax (by 
radiation) toward go = 0. The intercommutation con- 
verts these to NO strings! Figures [T] and |2| show frames 
from simulations of a right angle collision at V = 0.9, 
between go = 1 and go = 3 semilocal strings showing this 
behavior. We also considered collisions between unsta- 
ble (P > 1) Skyrmions. Again, intercommutation always 



occurred, and the reversion to NO always occured. The 
full range of the collision simulations is available online 
at http : //gravity . psu . edu/numrel/strings/ 

We have numerically studied collisions of strings and 
Skyrmions in the semilocal model. Our results demon- 
strate that collisions of any two objects in this model, 
whether they are identical or not, whether they are sta- 
ble or unstable configurations, leads to intcrcommuting 
for any set of collisional parameters. Additionally, the 
simulations all show a triggered collapse of the interact- 
ing objects toward a NO configuration. During the col- 
lision we also observe a large amount of energy lost in 
radiation. This may provide an additional signature for 
cosmological semilocal strings. 
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